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1.  INTRODUCTION 


The  use  of  digital  computers  to  solve  partial  differential  equations 
in  two  or  more  dimensions,  e.g.,in  fluid  dynamics  or  structural 
mechanics,  typically  involves  the  generation  of  numerical  values  for  the 
dependent  variable  (or  variables)  at  descrete  points  on  a surface,  which 
may  or  may  not  be  a plane.  If  more  than  a few  points  are  involved, 
graphical  presentation  of  the  results  is  virtually  mandatory  for  human 
comprehension.  Contour  plots  of  the  dependent  variable  convey  infor- 
mation very  effectively  and  hence  are  popular.  Two  algorithms  which  can 
be  used  when  the  descrete  solution  points  are  arranged  on  a rectangular 
grid  are  described  in  Reference  1. 

However,  not  all  problems  are  amenable  to  description  on  a rectangu- 
lar grid.  Finite  element  methods  frequently  use  triangular  elements  and 
some  solution  domains  simply  are  not  rectangular.  One  also  occasionally 
encounters  data  which  is  scattered  irregularly  over  the  surface.  The 
algorithm  described  herein  works  very  well  in  such  environments.  Though 
it  requires  a great  deal  more  memory  than  LEVEL  1 or  LEVEL  2 and  requires 
more  effort  to  set  up,  it  is  not  difficult  to  use  and  produces  high 
quality  plots  with  little  expenditure  of  computer  time. 

2.  THE  ALGORITHM 

2.1  Overview 

The  algorithm  presented  here  as  IS0VAR2  and  subsidiary  routines 
IS0VAR3,  FOLLOW,  DRAWBND,  DRAWL,  and  LABEL  is  a FORTRAN  recoding  of  the 
Algol  procedure  tricont  2 described  in  Reference  2.  It  retains  all  the 
capabilities  and  limitations  of  the  original.  Because  of  the  differences 
between  FORTRAN  and  Algol,  the  method  of  invoking  the  algorithm  is 
greatly  different  and  is  described  in  detail  below  and  with  some  examples. 

2.2  Operation 

IS0VAR2  requires  the  user  to  supply  several  arrays  and  simple 
integers  which  describe  the  arrangement  of  the  grid  and  the  data  on  it. 
Some  must  be  preset  by  the  user  while  others  may  be  treated  as  working 
storage.  We  shall  describe  them  in  order: 
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EL 

EP 

E 

PX,  PY,  PZ 


Integer  working  array  dimensioned  (E,  3). 

Integer  array  dimensioned  (E,  3)  in  which  EP  (R,S), 
(S  = 1,  2,  3,)  are  the  subscripts  in  PX  and  PT  of 
the  vertices  of  element  R,  proceeding  clockwise 
around  the  element.  See  Figure  1. 

The  number  of  triangular  elements  in  the  grid. 
Integer . 

Real  arrays  dimensioned  (G)  containing  the  data  to 
be  contoured,  X,  Y,  and  Z values  respectively. 

They  may  be  doubly  dimensioned  in  the  calling  pro- 
gram, in  which  case  G is  the  product  of  the  two 
dimensions.  Note:  PX  and  PY  must  be  scaled  to 
represent  actual  inches  of  pen  movement  on  the 
plot.  See  the  examples. 


The  number  of  points  in  each  of  PX,  PY,  and  PZ. 
Integer. 


HTS 


Real  array  containing  the  contour  heights  to  be 
drawn. 


K The  number  of  contours  in  HTS.  Integer. 

LL,  LP,  LR  Integer  working  arrays,  dimensioned  (N,  2). 

UNUSED  Logical  working  array,  dimensioned  (N) . 

N Maximum  number  of  lines  bounding  the  triangles. 

Lines  common  to  two  triangles  are  counted  once.  N 
is  never  more  than  3*G,  nor  less  than  E + G - 1. 
See  Figure  3. 
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Figure  1 

Definition  of  EP 

IS0VAR2  will  print  the  contents  of  the  various  arrays  if  the  logical 
variable  DBUG  in  labelled  COMMON  /BUGBUG/  is  set  .TRUE.  This  option 
should  be  used  with  discretion.  Large  meshes  will  yield  very  large 
quantities  of  printout. 

2.3  A Simple  Example 

Consider  an  array  of  30  Z values  arranged  on  a 6 x 5 grid  as 
shown  in  Figure  2.  X values  run  from  0 to  5 and  Y values  from  0 to  A. 

Let  us  create  a triangular  mesh  by  running  diagonals  from  upper  left  to 
lower  right  of  each  mesh  rectangle  as  shown  in  Figure  3.  Understand  that 
this  choice  is  completely  arbitrary  for  this  example  and  the  other  dia- 
gonal would  have  served  just  as  well.  In  other  applications,  circum- 
stance' may  favor  a particular  choice  or  combination  of  choices.  In  this 
example  we  could  have  taken  some  diagonals  one  way  and  some  the  other. 
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Figure  2 
Sample  Z Array 


We  can  then  create  the  PX,  PY,  and  PZ  arrays.  From  the  description 
of  the  problem,  G is  30.  Thus: 

DIMENSION  PX(6,  5),  PY(6,  5),  PZ(6,  5) 

DATA  (PZ  (I,  1),  I = 1,  6)  / 2.0,  3.0,  6.0,  7.0,  4.0,  3.0  / 

DATA  (PZ  (I,  5),  1=1,  6)  etc.  DO  1 J = 1,  5 

J1  - J - 1 

DO  1 I*  1,  6 

PX  (I,  J)  = I - 1 

PY  (I,  J)  = J1 

CONTINUE 

Then  the  EP  array  is  set  up.  From  Figure  3,  we  know  there  are  40 
triangles.  Therefore  E is  40.  For  our  own  reference  (The  information 
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is  not  passed  explicitly  to  1S0VAR2)  we  need  to  adopt  a numbering  con- 
vention for  the  triangles  and  their  vertices.  Let  us  therefore  start  in 
the  lower  left  corne  of  the  mesh  and  proceed  to  the  right,  numbering 
the  triangles  in  sequence,  with  odd  numbers  going  to  the  lower  triangles 
and  even  numbers  to  the  uppers,  as  shown  in  Figure  3. 


Cl.n  (2.1)  (3.1)  (4.1)  (5.1)  (6.1) 

Figure  3 


Triangular  Mesh 
Numbers  in  parentheses  are 
FORTRAN  subscripts  for  PX, 

PY,  PZ 

Numbers  in  triangles  count  the 
mesh  elements  and  correspond 
to  the  first  subscript  R of 
EP  (R,  S). 

We  also  need  a system  for  identifying  the  vertices  of  the  triangles. 
Therefore,  let  us  start  at  the  lower  left  vertex  of  each  lower  triangle 
and  at  the  bottom  vertex  of  each  upper  triangle  and  count  clockwise  around 
each.  c 
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All  that  remains  is  to  convert  the  double  FORTRAN  subscripts  to 
single  integers  for  EP.  For  this  we  use  the  same  method  FORTRAN  uses: 

The  relative  location  of  Azj  in  an  array  dimensioned  A(I,  J)  is  given  by 

L = i - 1 + I * (j  - 1)  (1) 

so  that  the  location  of  the  first  element  relative  to  itself  is  zero. 

Since  we  want  equivalent  subscript,  we  use 

L 1 = i + I * (j  - 1)  (2) 

The  entries  in  EP  are  the  equivalent  subscripts  of  the  three  vertices 
of  the  Rth  triangle.  Thus  EP(1,  1)  points  to  the  first  corner  of  the 
first  triangle,  EP(1,  2)  points  to  its  second  corner,  and  so  on.  In  cur 
example, 

EP  (1,  1)  = 1 + 6 * (1  - 1)  = 1 

EP (1 , 2)  + 1 + 6 * (2  - 1)  = 7 

EP (1 , 3)  = 2 + 6 * (1  - 1)  = 2 

EP (40,  1)  = 6 + 6 * (4  - 1)  = 24 

EP  ( 4 0 , 2)  = 5 + 6 * (5  - 1)  = 29 

EP  ( 4 0 , 3)  = 6 + 6 * (5  - 1)  = 30 

The  choice  of  contours  depends  on  the  problem.  We  will  use  the 

sequence  2.5  to  8.5  in  steps  of  1.0,  which  makes  K = 7. 

Lastly,  we  must  establish  the  size  of  the  arrays  LL,  LP,  LR,  and 
UNUSED,  i.e.,  the  value  of  N.  We  know  that  N is  not  more  than  90(=  3 * G) 

nor  less  than  69(=E  + G - 1).  If  we  count  the  lines  in  Figure  3,  we  will 

find  69.  For  simple  geometries  like  this,  the  lower  limit  will  always 
suffice . 

With  this,  we  are  ready  to  plot.  The  complete  driver  program  is 
listed  with  IS0VAR2  in  the  Appendix  and  the  completed  plot  is  shown  in 
Figure  4.  Note  that  IS0VAR2  draws  the  grid  boundaries;  no  special 
provision  for  those  is  necessary. 
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Figure  4 

Sample  Co.'.tour  Plot 

2.4  Multiple  Data  Sets  on  the  Same  Mesh 

If  one  wishes  to  plot  contours  of  more  than  one  variable 
over  the  same  mesh,  he  can  bypass  1S0VAR2  on  the  second  and  subsequent 
calls  and  enter  the  algorittmi  at  IS0VAR3.  The  arguments  are  EL,  EP,  E, 

PX,  PY , PZ,  G,  HTS,  K,  LL,  LP , LR,  UNUSED,  and  N.  Axl  have  the  same 
meanings  as  for  IS0VAR2  except  that  PZ  would  be  the  alternate  dependent 
variable  set.  The  user  must  establish  a new  plot  origin  or  the  new 
plot  will  be  superimposed  on  the  old.  Figure  5 shows  a sample  obtained 
by  "relaxing"  the  first  data  set  one  cycle,  i.e.,  for  every  point  not  on 
a boundary 

ZP  (I,  J)  = (PZ(I,  J + 1)  + PZ(I,  J - 1)  + PZ(I  + 1,  J)  + PZ(I  -1,  J))  / 4.0 
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Figure  5 

Contours  of  Second  Data  Set 

2.5  An  Example  from  Fluid  Mechanics 

The  primary  reason  for  developing  IS0VAR2  from  tricont  2 was  to 
plot  aerodynamic  data  at  the  exit  of  a turbine  blade  row.  The  general 
grid  arrangement  is  similar  to  Figure  3,  but  with  three  differences. 

The  upper  and  lower  boundaries  are  circular  arcs  and  the  left  and  right 
boundaries  are  identical  curves  separated  by  the  blade  pitch  angle.  The 
grid  is  typically  around  30  x 30  with  density  increasing  toward  each  of 
the  four  walls.  Figure  6 shows  a typical  plot  of  kinetic  energy  distri- 
bution, the  highest  levels  being  near  the  center. 

Figure  6 required  only  five  seconds  of  CDC  6600  central  processor 
time  to  generate,  but  the  total  array  storage  (including  PX,  PY,  and  PZ) 
was  more  than  31000  (decimal)  words. 

2.6  Potential  Improvements 
2.6.1  Memory  Requirements 

As  seen  in  the  last  example,  one  of  the  obstacles  to  the 
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Figure  6 

Kinetic  Energy  Contours  at  Stator  Exit  Plane 

general  application  of  IS0VAR2  is  the  large  amount  of  memory  it  requires 
for  the  working  arrays  EL,  EP,  etc.  Only  the  PX,  PY,  and  PZ  arrays  are 
floating  point  and  use  an  entire  computer  word  for  each  element.  The 
remaining  arrays  are  integer  or  logical  and  therefore  use  only  a fraction 
of  a word  per  element.  The  EL  and  EP  arrays  for  instance  are  dimensioned 
(E,  3).  These  could  be  packed  by  threes  into  a 60  bit  CDC  word.  Allowing 
18  bits  each  would  permit  a maximum  value  of  777777  octal  or  262143 
decimal,  large  enough  for  any  application.  The  LL,  LP,  and  LR  arrays  are 
dimensioned  (N,  2).  These  could  be  packed  two  to  a word  or,  by  storing 
corresponding  elements  of  LL,  LP,  and  LR  together,  three  to  a word.  Thus, 
half  to  two  thirds  of  the  integer  array  storage  could  be  saved.  The 
logical  array  UNTTSED  uses  only  one  bit  per  word.  These  could  be  bunched 
60  to  a CDC  word  to  save  more  than  98%  of  the  required  storage. 

Packing  and  unpacking  is  a simple  operation  in  CDC  Fortran  because  it 
includes  shift  and  mask  instructions.  IBM  360/370  Fortran  allows  a 
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length  specification  in  type  statements,  e.g.,  INTEGER  * 2 E,  which 
specifies  the  number  of  bytes  each  variable  is  to  occupy.  Other  com- 
puters may  require  assembly  language  subroutines  for  packing  and 
unpacking.  In  any  case,  the  modifications  would  require  only  a few 
hours'  work,  save  much  memory,  and  increase  execution  time  by  only  a 
very  small  amount. 

2.6.1  Labeling  the  Contours 

Many  times  the  plotted  contours  fall  close  together, 
double  back,  twist  and  turn,  and  generally  become  hard  to  distinguish. 
Figure  6 is  a good  example.  For  these  cases,  subroutine  DRAWL  can  be 
modified  to  mark  the  beginning,  end,  and  selected  intermediate  points  of 
a contour  with  a distinctive  symbol.  The  Calcomp  routine  SYMBOL  has  a 
repertoire  of  15  suitable  symbols  which  could  be  invoked  by  passing  the 
value  of  the  contour  loop  counter  J from  IS0VAR3  to  DRAWL  through  a 
labelled  COMMON  block.  Alternatively,  one  could  cause  DRAWL  to  print 
the  X and  Y coordinates  at  selected  intervals  on  the  contour.  The 
latter  method  is  used  in  the  program  which  generated  Figure  6.  Annota- 
ting the  contours  directly  is  possible  but  would  be  difficult  for  the 
very  reasons  stated  above. 

3.  CONCLUSIONS 

IS0VAR2  has  proven  to  be  a useful  and  versatile  algorithm  which 
produces  very  satisfactory  contour  plots.  It  requires  a lot  of  memory 
but  very  little  execution  time.  The  memory  requirement  can  be  reduced, 
as  described  in  Section  2.6.1. 
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APPENDIX 
PROGRAM  LISTINGS 

The  complete  FORTRAN  code  of  IS0VAR2  and  the  driver  program  TWO 
which  was  used  to  generate  Figures  4 and  5 are  listed  herein. 

The  listings  contain  occasional  continuation  lines  marked  with  a 
dollar  sign.  These  were  created  by  the  listing  program  to  maintain  the 
right-hand  margin  and  do  not  appear  in  the  actual  code. 
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PFOGP.AH  TWC  (FltT,  OUTPUT=4029,  T B r r 6=0uTr  UT) 
CI^NSION  w(7),  Ll(69,2>,  Lf(69,2»,  LM69,2>,  X(6,5), 
7 Y (6,  6)  , 

1 Z (6 , 5 ) f Z2  <6  » c ) 

LOGICAL  n°UGf  1(6°) 

INTCGGC  fr(«4C,Z),  FI  (40,5),  r,  C,  Q 
COMMON  /6UGpUG/  GBUG 


CAT« 

H / 2.5,  7.5, 

4. 

5,  5, 

5,  6. 

c,  7. 

5,  9. 

5 / 

DATA 

(Z (1,1 ) ,1=1,6) 

/ 

2.  0, 

3.  0 , 

6. 

7.0, 

4.3, 

r 

/ 

DAT* 

(Z(T, 2) ,1=1,6) 

/ 

4.  P, 

6.0, 

Q.% 

7.% 

6.0  , 

4.  C 

/ 

CATA 

(Zll,’)  » T=  3 > 6) 

/ 

4.  ", 

6.0  , 

P.r, 

5.0, 

6.0  , 

5.  0 

/ 

DATA 

(Z(I,4) ,1=3,6) 

/ 

3.0, 

6.  C , 

7.C, 

6.0, 

6.0  , 

3.  r 

/ 

DATA 

(Z(l, 5) ,T=1,6> 

/ 

3.0, 

4.0, 

6.0, 

5.0, 

4.0  , 

2.  C 

/ 

K ( I , J ) 

nr  1 

FI 

DO 


= I ♦ 6 * (J  - 


DO  i j 

X (I,  J) 

Y (I, J) 

contikuf 


OPUG 

F 

6 

N 

0 


1,  6 
3 - 1 
= 1,  5 
= FI 
= J - 1 


= .falsf. 

= 4 C 
= 3G 
= 6C 
= 7 


L 0 AO  pc  4JTH  X,  Y,  AWT  7 SUP^CPIPTS 


CO  ? I 
LiO  ? J 


= 1,  5 

= 1*  4 


L 0 W r 4 TFIAKGLE  HAS  SUBSCRIPT  F TN  FP  : UF=»FR  HAS 
SUeSCFIFT  L. 


FO  <H,2> 
EP  (M,3) 

E o <1  » 1 > 
EP  (L  ,2) 
FOCL,3) 
FONT INUF 


=2*  (I  ♦ 5*  (J 
= L - 1 
= K(T,J> 

= <<I,J  ♦ 13 
= K<I  ♦ 1, J) 

= K ( I + 1,J) 

= K(I,J  ♦ 1) 

= K ( I ♦ 1, J ♦ 1) 


FALL  PLOT  (2.0,  t.0,  -3) 


SB 

I* 

z S 

in  c/> 


I 


CALL  ISO  VAR2  (tl,  FP,  f,  x,  v,  z,  r. , h,  Q,  LI,  LP, 
* LR,  U,  N') 

C 

r GPNFPATt  SFf  CNO  TATA  i>FT 

C 

no  3 1 = 1,  6 

Z?(I,1>  = Z(T,1) 

7?  ( T » 5 ) = Z(I,S) 

3 continue 

no  4 J =7,4 

Z?(1,J)  = Z(1,J> 

Zu  < ft  * J ) = 7(6, J) 

A fONTINMF 

nr>  5 T = ?,  5 

JO  5 J = 2,  u 

= l.Zc  * (Z(I,J  «•!>«■  Z ( I , J - 1)  ♦ 

« Z(I  1,J)  ♦ 

1 Z (T  - 1, J) ) 

5 n ON-r  JF'UE 
C 


CALI 

°LCT  (10.", 

O.C,  -3) 

CALI 

xsrh/^pj  (x, 

Y,  Z2,  r„ 

0,  U , LF,  LR,  U, 

* ) 

CALL 

PLO”  (P.5, 

f'+.C,  -3' 

CALL 

S /‘"FPL  (r.r 

, ".5,  0 . 1 n F , 

•*cIK'ISHF.O“,  Sr.  0, 

h ) 

CALL 

PLTTt 

ST  OF 

r«jr 


SUn*DU7TN£  IFCwAP2  (EL,  EF,  E , PX,  Fy,  PZ,  G , HT  F , 
* K,  LL,  IF, 

1 L9i  GNUSFE,  K) 

COMMON  /BUGdLG/  G»UG 
CCMMON  /ISOLIN/  L,  LP 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

r 

C 

C 

C 

C 

C 

/* 

c 

c 

c 

c 

c 

r 

C 

C 

r 

C 

c 

r 

C 

C 

C 


T H IF  A LGCF  77  HN.  IS  USrC  FO^  C ONT  CURING  OVER  AN 

* I»FEGULAF, 

TPIANGULAP  MESH.  Thp  P^rAHTTFrc  FFOUlFtFO  ARP  ! 

F - THE  NUN6FP  OF  TFIANGULAF  rL  rM  t FT  S • TNT cGtP. 

El  - INTGEF  WORKING  AFfAY,  C lKr,''STONPrj  <E  , 

FP  - INTEGER  AF  F A y 01  Hr  K’SI  0»‘f  0 (r,  7),  IN  N'HICF 
« FC(P,  S), 

(S  = 1,  2,  7),  mPl  THF  SUE'* OF  1PT S IN  FX  ARC 
« P*  CF  TR£ 

PCPNERS  OF  THE  ElcNEN”  n,  APR ANGrG  CLOCKWISE 
AfCLNC  THF  E'l  f HE  NT  , 

PX  - kFAl  APfAYF  CONTAINING  THf  "ATA  TO  5E 
$ CCNTOLFcr# 

OY 

<57 

G - T Hc  NUWErc  CP  PflN’S  IN  EACH  OF  PX  , PY,  PZ. 

* TNTEGFC. 

hTS  - peal  AFFAY  of  CCn'tOUP  HETCHtS  TP  Bt  ORAWN. 

K - the  NUMFEF  CF  CO^TCL'PS  IN  H-S.  integer. 

LI  - INTPGFF  WORKING  ARRAYS,  DIMENSIONED  (N,  2). 

LP 

IP 

UHUSEn  - LOGICAL  WCiPKIN’G  AR^AY,  cTMEnsTuNF")  (N)  . 

N - MAXIMUM  KUHpfR  Cp  LINT  «C'JNDING  THE 

* TRIANGLES. 

VFVFK  corf  than  ? * G,  nor  UPSS  THAN  E ♦ C - 
? 1. 

F XA»'P ir 

FOR  A.  PFtTA  N’GUl  A n AFT-Ay  CO  NT  AIDING  rt  PLEHcMS  IN 
f T HE 

X OIRfOTICN  ANC  0 EL  rM£NTS  IN  ThF  Y DIRECTION, 

F = 2 * <h  - i>  * CO  - 1) 

G = ► ■»  0 
K = YPUC  CHC.ICP 
N = F * G 

A PFTAUFO  CFSCPIFTTON  CF  THF  ALGORITHM  IS  GI\,cN 
$ IN 

SECTION  3 OF  THE  ORIGINAL  PEPpo". 

AOAFTPT  FRC*  TmF  ALGOL  frnr^CUrr  TPTOONT  ?,  FCUND 

* iv 
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ooc  oi  tansiwu  — * «, 

ngronoraa  uimt  jsjs  rr  tsyj  sihj 


rwiiinond  i r»  o .*»  r»  r»  o ■“>*-»  o ->  ->  o i->  -»  <-> 


( iP  1-K1C-1D  ALPOcI^HHS  TMr  p POOL)''  T 10  K CF 

T roKTOttF  hafs, 

9.W.  Hr  AF  <NATl.  PHY**.  1*9.,  7 rn0TNGT0N , c«r»L3»ir) 

uAPcM,  ly??.  PAG*  21. 

TnHT/LLY  7 HF  PulINnfiry  LT,T  '•C"kjT,  L3,  ANC  THE 

* PLP**FK7« 

0 r 1 Ht  AFP  AY  El  APE  S-T  T0  ZF^C  AFC  THC  T f) T A L 

* LXNF.  F CUNT  , 

L,  If  SET  Tc  JT  MAXTHIIF  pn?ST'3LE  VALUF. 


DIKE  N^ICN 

HTS  <K), 

LL 

(N,2)  , LF 

<N»2>* 

LP  (N,2), 

s 

py  (G),  fv 

1 

<G>,  p 7 

( G > 

I Nr  LGFu. 

E,  El  (E 

» 3 ) , 

FF  ( F , 1 ) , 

G*  PA, 

PB  , PC , F , 

? 

S,  FA,  PP 

LPG1CAI 

CPL'G,  UN 

USE  r 

(Ml 

Lc  = 3 * F 

L = LF 

19  =0 

pP  1 «?  = 1,  F 

lIO  1 F =1,3 
rl  (P,  S)  = 0 

1 7QNTJ  NUF 

rACM  FI  F HE  NT  OF  THE  *»E  S H IS  EXAPINEO  IN  TUAN,  TTS 

* foist 

LAdtLS  EXTEACTFu,  Ant  ’HP  AFPECFFTATF  ENTFJFS  tk 

$ TU|C 

aivt^Y  el  MAfF  HSTNG  SUFROUTINF  LA  EEL. 

nn  2 R = 1,  F 

p*  = fc(F,  u 

pp  = F r < p , 2) 

rC  = EP(H,  3) 

iF  < EL  <f  * 3 > .FQ.  C ) TALL  LABEL  (r,i,P3,pC,F,PF,EL> 

IF  <*l<«?,2>  .ED.  0)  C * LL  LAQH.  ( p , ? ,=>C , P A , F , F P , E L ) 
IF  (FL(i%,3)  .EO.  A)  CALL  LA9FL  < » , 3 , p A , p E,  F , E°  , E L ) 

2 roNTmip 


THP  \ Inf  CCL'VT,  L,  1«  AVFNOPF  A K'L  THP  L IN  t LA^IS 
flnju^Tcj  so  THA”  T H E SOUNTAPv  tTNES  HAVE  LA  EL' 

? It  2*  3 , 

...,  Id,  ANC  the  TNTFCICP  LTNE4,  have  LABELS  L?M, 
f L p *2 , 

• • • 9 L • 


n o n r>  o n o o n n n n o 


» 


►'* 


L = l - LC 

no  3 p = j , F 

00  * S = 1,3 

I p Cl  <P,S>  .or  . IR)  Fl  <F,<n  = EL  <P,S)  - 1C 
3 nONTTNLE 

TH^  ni^FNSDCNS  OF  Thr  AFFAYS  IP,LI,L»  <L,2> 

♦ FFOUIFFC 

FOP  ISPVAF?  «Ft  NOW  KNOWN. 

It'T’TALLY  ALL  r N T F 1 1 c T ► T Hf  AREA  VS  Lh  AND  LL 
A OF  SF1  "0  ZrtO. 

h CONTINUE 

IF  (u  ,LF  . N ) C-C  to  5 
PRINT  " , "TNCRFASF  N TO  ",  L 
5- TOP 

5 CCNTTNUf 

DO  6 r = 1,  L 

DO  f c = 1,  2 

l.R(F , S)  = LL  <P,  S)  = n 

6 OONTINUF 


THf  ARFAYS  IF,  Lr,  AND  LL  APF  FFT  UF  PY  SCANNING 
$ THRCLCH 

ALL  FLFMFNTS  OF  THF  MF^H  ANO  FXTPACTINS  THE 
$ A FFROFF IAT  E 

FNTPIFO  CF  ARFAYS  EF  ANC  FL. 


7 


b 


DO  7 p 

= 1.  F 

(/ 

o 

c 

= 1,  3 

M® 

= LI  (n, 

D 

SA 

= * + 1 

IF  (S 

.FC.  3)  SA 

= 1 

SO 

- S - 1 

IF  <S 

.EO.  1)  Sf? 

- 3 

NA 

= 2 

IF  (LR (M*  , 1)  .LG. 

0)  NA 

LRIWA, 

NA ) = FLIP, 

SS) 

IL(WA  , 

NA)  = El  IP, 

SA) 

IF  (NA 

.n*-,  i)  GO 

TO  7 

I P (HA  , 

1)  = EP(P, 

CA> 

LP  (MA  , 

2)  = F P < R , 

S 0 ) 

CONTINUE 

IF  ( • NC‘r , 

C0UC-)  GO  TO 

12 

phint  ion 

no  n i 

= 1.  c 

PPINT  101,  I,  (EP(I 

, J)  , 

rONTTNUr 

PRINT  1C2, 

L,  l.« 

ppINt 


in; 


1 


J = l, 


T) 
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o o o 


O'  IN’  lFA,  F,  IF  (Fi,  1),  PX  ( Lc  ( R » 1)),  P Y (LF 

$ (P,  1)  ) , LF 

1 (F,  2),  FX  <t_F  (r,  2»),  rv  ((  r rr,  2)) 

Q CONTINUE 

PxJN’  3 F 5 

PO  ir  r s 1,  L 

P=>TNT  101,  F,  LL  (P,  1),  U <F,  2),  Lc  R,  1),  LR 
S (P,  2) 

13  CPnTTnmf 

Pfc'INT  10  E 

no  li  c = i,  f 

PPINT  !U,  F , <n  (P>  S>,  c = l,  ■»' 

11  f'ONTINUF 

12  CCNiTNUF 

u ONT  OURF  Aft  NOW  CPA  WN  USTKF  SUcrCUlINt  ISCVAP?. 

CALL  T<?°VA  F3  (PX,  PY  , F7,  G,  H7S,  * , LL,  LF,  LR, 

♦ UNUS  c 0 , N) 

R FT  U ° N 


100 

FCpM  at 

< 

1 HI  , 

AX, 

1 HI , 3 X , 

7HFF  (1,1), 

3X,  7HPF(],2) 

$ 

, 3X  , 

1 

7HFF  ( 

3)  3 

101 

FORMAT 

( 

1H  , 

5iir 

) 

102 

FORMAT 

( 

*0 

L =* 

, I?,  * 

Lfi 

=*,  13  ) 

103 

FPKM  «t 

( 

*1 

LTMF 

LIST*  / 

7X, 

AHl.IXe, 

bX,  1.MPO  CF, 

? 9 X,  1 H X , 9X, 

1 1HY,  6X,  2HTO,  9X,  1HX,  9X,  1 H Y ) 

FORMAT  ( ih  , 110,  2(110,  ?F10.1)  ) 

ppruftT  ( ♦ l LI Nr  LI^<SGFF*  / 7V,  LHLINF,  3X, 

i 7HLL  (P,  1)  , 3X , 

1 7Hl»  (R,2),  7X,  7HLf  (SI),  7*J<h(R,2)  ) 

FOFMA  t ( *1  CCONCAPY  LINES’  / 7X,  8HTi-  TANGLE,  3X, 

i 7HFL  (°, 1)  , 

1 3X  , 7 HFL (F , 2 ) » 3X , 7HrL(R,3)  ) 


loonnnonn  n n .*)  o n o n n n o o n n n Tnoonnnnon 


i 


bUBROUTTKE  TSCVAf?  (PX,  FT,  F7,  G,  HTS,  K,  LL,  LF, 

? LR,  JNUSFC, 

1 N) 

CCHMON  /OUGBUG/  n»ur, 

common  /flips/  open,  first,  last,  h 

C CMrl  OKI  /IFtLTH/  L,  Li 5 

THIS  FTCCF r LFE  JS  A SUBFf-oerDU^E  0 F IS0VAR2  9LT 
« raw  BF 

USED  SEPARATELY.  IT  IS  A L SC  FOP  CONTOURING  OVER 

<t  f m 

IPPctULAF  TF IANGULAR  HESH  • THr  PAr AHFTERS 
$ PFQUIKFr  APE 

Lit,  T He  NURBEF  OF  90UNCAFY  LTNPS  IN  THE  MESH, 

L,  the  TOTAL  NUHBEF  Oc  LINES  IN  THE  MESH  , 
f TNTfcr,rf  4 f . r> A Y S 

LP,  LP,  IL,  LhERr  L c ( " , N ) , (N=l,2),  ARE  THE 

* L flcFLS  CF  THE 

POINTS  AT  THE  Fpns  rp  tlNr  N , AND  THE  EMMIES  IN 

* L;?(P,N> 

AND  LL(H,N)  SPFdFY  THE  L«BtLS  OF  THE  LINES 
? JOINING  A T 

THFi,6  POINTS.  ALSO  WEEC^P  ARE  " Hr  ARRAYS  FX,  CY, 
T t ME  C7 

AND  THE  INTEGER  K,  W^IO  ARr  AS  F OF  ISO  VAR  2, 

* FULL  CETATLF 

ARF  GIVEN  IN  S E C T 1 0 N 7 OF  THE  ORIGINAL  REPORT. 

T HE  ocUNDAFY  OF  THE  »*ESH  TS  NOW  DRAWN. 

DIPENFICN  HTS  (K>,  u ( N ,2 ) , Lf  (N,2),  LP  <N,2), 

* PX  ( G) , PY 

1 (G),  FZ  (G) 

INTEGFF  G,  FA,  FB,  TEMP 

LOGICAL  CcLG»  UNUSFD  (N) 

CALL  OPAWHNn  tFX,  PY , G,  LF,  N) 

EACH  l T Nfc  IS  E X A N 1NC  D IN  X(jrN,  THE  HEIGHTS  OF  ITS 

* ENC 

POINTS  extracted*  anc  thf  entptls  IN  the  ARRAYS 
S IF,  LS 

AND  LL  INTERCHANGED  IF  NFCFSSAJ'Y  SO  THAT  C0NTCUFS 
S A rE 

TRACED  WITH  HIGH  GROUND  ON  THF  RIGHT. 


no  i h 

« 1*  L 

PA 

= LF(m, 

1) 

P9 

= IF(M, 

2) 

ZA 

= F 7 < FA ) 
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o o o o n o o n nn  n o o n n n 


zq  = rz(HB) 

IF  ( 7 a .GT.  Ic)  GO  TO  1 


L° ( w , 1) 

= 

p .1 

in  (ii,  ?) 

= 

FA 

TEHF 

= 

L°  (*', 

1) 

LR(“,  1) 

= 

L F ( * , 

2) 

. LP  ( *'  , 2) 

= 

TENP 

»f  Mp 

rr 

L L ( M , 

1) 

LL  (M,  1 ) 

r 

LL  (*, 

2) 

l L ( * , 2 > 

r 

T Pf  P 

FONT  I N'Ur 
IF  ( .NO'. 

ceuco  go 

TO 

PPIN ' 101 

on  2 “ 

- 

1*  l 

PR  TNT  101,  E , LF  (“ 
? ( M , 2 ) , l P 

1 <-,,  1),  ur  (“,  2) 

2 PONTTNUE 
PRINT  1C  2 

3 CCNT  TNUr 
C 

0 FAOH  rCNTQUf  HFIGHT 

C 

no  6 j - i,  < 

H = HTS  f J) 


1)  , L r C , 2)  , LI.  ( 1 , 1)  , LL 


T S Fr  ALT  WJ'|  M IN  TUFN. 


THF  AE  r A Y IN  USE  0 IS  S>ZT  UD  FOr  '.H3S  HFIGHT, 
p L\»ussr(“) 

q e X N G SFT  . 1 «=(!tr  # jp  ^ PtNTOUF  1 TNE  Or  THIS  HEIGHT 
OROSSFS  LINf  P « 

00  4“  = 1 , L 

UNUSE  0 ( m)  = PZ(LP(M,  1)) , GE.H, ANC.PT ( L F ( “ , 2)) 

$ .LT.H 

4 PON'TINUf 

THF  AFFCy  UCSFO  IS  SGANN^r  F OP  THE  START  CF  A 
T CCNTOl'F 

AND  THE  CONTOUF  IS  tcaCEO  USINF  SUBROUTINE 
» FOLLOW. 

no  5 N = 1 , L 

WE  USE  I TKciroo  of  M MFPE  BECAUSE  FOLLCW  RESETS 
T H Tc  If'TFGFE  AND  WOLIC  “rSS  Uc  the  LOOP  .TCUNT. 


X = w 

if  (UNusrctM  .and.  l°  if , i)  , a' r « t)  oul  follow 

* (T,PX,PY,P7,G 
1 , LI  ,LF, I F, UNUSED, >') 
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101 

102 


S 7HLI  (M, 2) , ?X, 

2 7HLR(M,i),  2X,  7HLR(V,2) 


FORMAT 

FORMAT 


( lb  , 
( 1 HI 


71  in  ) 
) 


FNO 
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oooooo  o r->  o o o o.-iooo  oooooo 


? L9n  OUT  IN  F FClLCiW  <«,  ex,  PY,  rZ>  ~ , LI,  L°  , L5, 

$ LINUS'"  0,  m ) 

rc^MOM  /SUGEIC/  C ^ 1 1 G 

CCMUOS  /tccIJ^/  L,  LP 

POHON  / c L UF  c / OPEN  , r I F ?T  , L A c T , H 

this  FucftcinjN17  is  usfp  "o  follow  a particular 
* CONTOUR 

THROUGH  "r  H r VCFM.  THF  F AF  AN  rT  r " N 13  THE  LAPEL 

5 nr  T h f 

IT  Ml  O'  WHICH  CONGOU  ?EGJNC. 

DTMrNF7nK  LL  (N,2I,  L n (A,?),  IP  (U,2),  FX  <G),  F 
f (G>,  °Z  (0) 

T NT  E.  G c F G,  rf.  , P«,  TFvp 

LOGICAL  t,F  10,  cTi>ST,  l r?l,  orru,  UN'JSED  (M 

NSTAFT  - 0 

T HP  pr  COE  OLF  F BEGINS  ^ SrT'ING  THE  INITIAL 
« VALUrS  OF 

first,  laft,  ANT  OPEN, 

FtRFT  = * T r L F * 

l AFj  = .FALSE. 

OFFN  = K.Lc.lc 

ic  ( opug)  pf i nt  ina,H 

NEXT  F C T !•' T . 

r H c F C.  T N T WHFFF  ■rHF  lixr  f'FCcSF  S T Ht  L T N 7 V I c 
OFTF  r MTNFO  PY  TNYEPSr  T N r A F TN"  ER.PPLA  TIDN. 


1 fa 
P9 


= LF<i<,  II 
= Lr  ( * t 2) 

= F 2 ( c A ) 

= n/  (ro) 

= o.r 


ON 

= 7A  - 

ZP 

TF 

(ON  . 

ne.  r.in 

T(  — 

(H  - 

ze>  / on 

TP 

= l.r  - 

T A 

YF 

= Tf  * 

PX(PA) 

♦ Tc 

* CyfTC) 

YF 

= TA  * 

r Y ( n ft  > 

+ To 

* P Y ( P 3 ) 

IF 

( Q°UG) 

FFTN’ 

1CL , XP 

,YP 

T M F ftpoporrJflT 

t r 

L c NTP 

Y T N 

THr  A cr  A y UNUSEO 

. c a LS F 

V 

. UNI  LSS 

IT  is 

r he 

F Ir  ST  ^0Int  uF  A 

CONTOUF.  IF  IT  JC  f HF.  F T i ST  FO.TN1  OF  h OCNT0lR, 

? TMC 

L*BFL  CF  -Ff  INITIAL  L I Nr  Tc  ST0PFP.  OTHERWISE, 
<S  TrrTr 


k 


I 


apf  rAppirr  cut  to  scf  tf  t i*  the  last  pcint  cn 
the 

CONT  PUP. 

UNUSFP<*>  = FIRST  . AND.  .hOt.OdEN 
TF  (FIPST)  GO  TC  ? 

LAST  = (CFFN.AND.M.lE.l.e) 

.OR.  ( .NO""  .OPEK.ANO.M.FQ.mstact) 

GO  T>  3 

MST4PT  = m 

THE  COORDINATES  9F  T AFr  OUTPUT  TO  THtr  POOCEDLFE 
OF  AWL.  N 

CONTINUE 

CALL  PPAWl  <*F,  Y r ) 

DROVT  EE  D THAT  IT  IS  HOT  THE  LAST  POINT,  THE  L A EEL 
OF  t H F 

N - XT  LINE  THAI  THE  PONTpilo  CPOSSFS  IS  OETFCJhIAEC 
A N P THE 

POUT  THE  OOFS  PACK  TC  FIFO  A NEW  POINT  T. 

FIFST  = .FALSE. 

IF  (LAST)  RFTITN 

TEMP  = LF  (K,  1) 

I i caff  up  FKC  UP  Ik  rHL  BOCNIFP 

IF  (tfmp  .FC,  C)  PETUPN 

M = U (“,  1) 

IF  (UhUSEL'<TFNF)  ) K = ~EMF 
GO  T0  1 

F OF M A T ( *P  X A NO  V rO°  COMOt.,c*,  G15.7  ) 

FCPMAT  ( 5X,  EG  1 5, 7 ) 

END 


n o r>  o r>  o n o 


SUBRO'JTINF  PFA  K8Nr'  (PX,  RY,  G,  Lr  > N) 

C CUM  ON  /ISOLJP/  l,  L.P 

THI°  rnoCFCUPF  TC  USFO  TP  DP  t W THF  BOUNDARY  OF 

♦ TFF  mfsh. 

FACH  erUNDATY  LT*'L  IF  FCANNFD  1 N -URN,  THF 

* C CORCIMTFS  OF 

ITS  «MO  POINTS  FXTRAP’Br,  AN?  THF  A op  OPR  I ATE  FEN 
HOVFMFNTS  ►'ADE. 


1 


GTmENcIGN 
I NT  r GFF 
00  1 t 
"A 

CALL  PLOT 
PA 

F a ll  CLCT 
CONTINUE 
PFtUPN 
ENP 


LF  <N,2>,  FX  <G> 
G,  Fa 
= 1,  I p 
= L r ( I » 1) 
(FX(PA),  PYCPA), 
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